Bound states at the interface between antiferromagnets and superconductors 
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We present a detailed theoretical investigation of interfaces and junctions involving itinerant an- 
tiferromagnets. By solving the Bogoliubov-de Gennes equations with a tight-binding model on a 
square lattice, we study both the self-consistent order parameter fields proximate to interfaces be- 
tween antiferromagnets (AF) and s-wave (sSC) or d-wave (dSC) superconductors, the dispersion 
of quasiparticle subgap states at interfaces and interlayers, and the local density of states (LDOS) 
as a function of distance from the interface. In addition, we present the quasiclassical approach to 
interfaces and junctions involving itinerant antiferromagnets developed in an earlier paper. Ana- 
lytical results are in excellent agreement with what we obtain numerically. Strong efi'ects of pair 
breaking in the presence of low-energy interface Andreev states are found in particular for AF/sSC 
interfaces when interface potentials are not too high. Potential barriers induce additional extrema 
in the dispersive quasiparticle spectra with corresponding peaks in the LDOS. Discrete quasiparti- 
cle dispersive levels in AF - normal metal (N) - AF systems are found to strongly depend on the 
misorientation angle of the magnetizations in the two antiferromagnets. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



I. INTRODUCTION 



Interfaces of magnetic materials with normal metals 
and superconductors have attracted much attention in 
recent years because they can strongly influence prop- 
erties of mesoscopic and nanoscopic systems, and may 
play important role in compounds with competing mag- 
netic and superconducting ordering. Hybrid supercon- 
ducting systems involving ferro- and/or antiferromagnets 
manifest unusual properties associated with spin and or- 
bital effects, and are of both fundamental interest and 
important for technological applications. Ferromagnetic 
layers can spin polarize quasiparticle currents and Zee- 
man split surface densities of states, with possible ap- 
plications in spintronics^. Superconductor - ferromagnet 
-superconductor (SC/F/SC) junctions have been shown 
to display — it transitions with varying temperature, 
width, or orientational structure of the magnetization of 
the ferromagnetic intcrlayer^^-'^^. 

There are also many situations of fundamental and 
practical interest which involve interfaces with antifer- 
romagnets. In particular, many of the properties of 
high-temperature superconducting (HTS) cuprate mate- 
rials are thought to result from a competition between 
antiferromagnetic and superconducting order, and there 
are many naturally occurring situations and possible de- 
vices which might involve such boundaries. These in- 
clude interfaces of insulating and highly doped cuprates 
or superconductor - antiferromagnet - superconductor 
(SC/AF/SC) junctions"'i2^ HTS grain boundaries'^ 
where antiferromagnetism may play a role as a surface 
state, and the antiferromagnetism which has been ob- 
served in HTS vortex cores"'^"'"'^. At the same time, 
there exist only very preliminary results of experimental 
and theoretical investigations of proximity and Joseph- 
son effects through various types of antiferromagnetic 



interfaces'^"^^. Below we study theoretically interfaces 
between itinerant antiferromagnets and normal met- 
als or superconductors. Itinerant antiferromagnets like 
chromium and its alloys^°'^^'^^ are metals above the Neel 
temperature. In the antiferromagnetic phase, however, 
an energy gap in the quasiparticle spectrum arises ei- 
ther on the whole Fermi surface or on parts of it. Simi- 
lar properties are also manifested in Mott antiferromag- 
nets, in particular undoped cuprates. Since they possess 
strong correlations due to large on-site Coulomb repul- 
sion U , the mean- field approach of the present paper can- 
not be applied quantitatively to Mott systems, whereas 
it applies well to itinerant antiferromagnets with com- 
paratively small U. Wc expect, however, that our main 
conclusions regarding AF/S and AF/N interfaces can be 
qualitatively applied also to interfaces with Mott antifer- 
romagnets because they are based largely on symmetry 
properties and general characteristics of antiferromagnets 
like the doubling of period, nesting conditions and the 
wave vector of the antiferromagnetic pattern. In itiner- 
ant antiferromagnets, the energy gap in the quasiparticle 
spectrum is determined by the antiferromagnetic order 
parameter, i.e. the sublatticc electronic magnetization 
TO. This applies both to commensurate antiferromagnetic 
phases and to phases with spin-density waves^^, and is 
reminiscent of the situation in superconductors, where 
the energy gap is determined by the superconducting or- 
der parameter A. 

Recently, it has been demonstrated theoretically that 
a new spin-dependent channel of quasiparticle reflection, 
the so-called Q reflection, takes place at interfaces be- 
tween itinerant antiferromagnets and normal metals or 
superconductors^^. Parallel to the interface, the momen- 
tum component of low-energy normal-metal quasiparti- 
clcs changes by a wave-vector Qn in a Q reflection event, 
where Q is the wave- vector of the antiferromagnetic pat- 
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tern. Assuming small Fermi velocity mismatches and 
taking into account the nesting condition Ep{p + Q) = 
~Ep{p) in itinerant antiferromagnets, one can sec that 
a normal-metal quasiparticle changes its total momen- 
tum by Q and the respective velocity changes its sign 
in a Q reflection event. Hence, normal-metal quasipar- 
ticles with energies less than or comparable to the an- 
tiferromagnetic gap experience spin-dependent retrore- 
flection at antiferromagnet-normal metal (AF/N) trans- 
parent interfaces. Furthermore, Q reflection processes 
generate quasiparticle bound states below the AF gap in 
AF/N/AF junctions, analogously to the case of subgap 
states in SC/N/SC systems formed by Andreev reflec- 
tion. The AF/N/AF bound states arise from a coherent 
superposition of electrons with momenta k and k-l-Q and 
almost opposite velocities. Subgap states arise also at 
AF/SC interfaces as a combined effect of Andreev and Q 
reflections. Among a variety of subgap states, low-energy 
states with energies Eb ^ min{m. A} are of special in- 
terest since they can result in low-temperature anoma- 
lies in the Josephson critical current, as well as low-bias 
anomalies in the conductance. Low-energy quasiparticle 
interface states were also predicted to occur on antifcrro- 
magnetic - s-wave superconductor (AF/sSC) interfaces in 
the absence of specular reflection, when one can disregard 
effects of interface potential barriers and Fermi velocity 
mismatches. For an sSC/AF/sSC junction, these bound 
states are split due to a finite width of the AF interlaycr 
and carry the supercurrcnt. At AF/dSC interfaces, low- 
energy bound states Eb ^ min{m, A} do not exist, at 
least if the order parameters are small compared with the 
hopping matrix element (A,m <^ t). This is contrary to 
the case of a (110) surface of a dSC confined with impen- 
etrable potential wall where zero-energy surface Andreev 
states are formed^^"'^"'^ . 

Below, we extend the study of effects of Q reflec- 
tion processes based on self-consistent solutions of the 
Bogohubov-de Gennes (BdG) equations. This goes be- 
yond the framework of the preceding paper^^ . In general, 
we find excellent agreement with the results by Bobkova 
et aP'^ . At the same time, the more general approach 
of the present paper allows us to study several impor- 
tant new problems. In particular, we discuss the effects 
of interface potentials on the dispersive quasiparticle in- 
terface states and the LDOS. In the presence of poten- 
tial barriers and Fermi velocity mismatches, there exists 
an interplay of specular and Q refiections. We demon- 
strate that potential barriers on interfaces between AF 
and either s-wave or d-wave superconductors can result 
in new cxtrema in dispersive quasiparticle spectra and 
additional associated peaks in the LDOS. We also find 
new interface quasiparticle states with subgap energies 
near the edge of the continuum, which arise due to self- 
consistent suppression of the order parameters near the 
interface. By studying effects of interface pair breaking at 
AF/SC interfaces, we find that for the (110) orientation 
the self-consistent suppression of both antiferromagnctic 
and superconducting order parameters near the interface 



is accompanied by even-odd spatial oscillations. We show 
that discrete quasiparticle dispersive levels in AF-N-AF 
systems strongly depend on the relative orientation of 
the magnetizations in the two antiferromagnets. Effects 
of the misorientation angle turn out to be analogous to 
the influence of the phase difference on the discrete quasi- 
particle spectrum in SC-N-SC systems. 

The paper is organized as follows. In Section II, we in- 
troduce the microscopic model used to study various in- 
terfaces with antiferromagnets, the BdG equations with 
a mean-fleld treatment of both magnetism and super- 
conductivity general enough to study both s- and d-wave 
pairing symmetry, as well as various interface potentials. 
In Section III, we sketch the derivation of the associated 
quasiclassical (Andreev) equations, complemented with 
boundary conditions, by assuming slow spatial variations 
of both order parameters. In Section IV, we study quasi- 
particle states at (100) and (110) interfaces of both s- 
and d-wave superconductors with antiferromagnets, and 
compare the results of numerical evaluations of the BdG 
equations with the predictions of the quasiclassical the- 
ory. We end section IV by showing how the bound state 
energies can be obtained within a transfer matrix formal- 
ism. In Section V, we study the AF/N/AF junction, and 
discuss a novel "spin-7r" configuration where the relative 
phase of the staggered magnetization on both sides of 
the junction can tune the energy of the interface bound 
states. Conclusions and perspectives for future work are 
presented in Section VI. 

II. MODEL 

For studying the electronic structure of interfaces be- 
tween antiferromagnets and superconductors or normal 
metals, we consider the following 2D mean-field Hamil- 
tonian on a square lattice 

- ^{^i - hi)nicr + ^mi{hi^ - fiii) . (1) 

i(7 i 

Here, Ay and rrii denote the superconducting and mag- 
netic order parameters, respectively, c^^^ creates an elec- 
tron of spin a on site i, t denotes the nearest neighbor 
hopping integral, /x is the filling factor, and hia = cl^Cia 
is the particle number operator on site i. In Eq.(l), hi 
is an interface potential. We will study self-consistently 
only singlet s-wave or d-wave superconducting pairings 
defined as A^ = —{Vi/2)(ciiCj^ — Ci^cji). For s-wave 
pairing, one should put i = j, whereas the d-wave or- 
der parameter Ay connects nearest neighbor sites. The 
self-consistent magnetic order parameter is represented 
as rrii = {Ui/2) [{fii^) — (nil)]- In the bulk of the an- 
tiferromagnet and in the absence of any perturbations, 
the staggered magnetic order parameter takes the form 
ruj = (— l)-'°+^''m = cyip(iQj)m. For a square lattice 
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with the crystal coordinate axes a and h, the antiferro- 
magnetic wave vector is Q = (7r/a,7r/a). Within the 
framework of a generic Hubbard-hke model, the stag- 
gered antiferromagnetic gapped state is stable only at or 
near half filling. For this reason we assume below van- 
ishing or small ^. 

We choose a coordinate system where x and y describe 
coordinates perpendicular and parallel to the interface, 
respectively. For a (100) interface the x and y axes coin- 
cide with the crystal axes a and h. Then the normal-state 
electron band is given by 



^(k) = — 2t(cos ka 4- cos kb) — /i. 



(2) 



and the respective Brillouin zone is spanned by ka^h G 
[— TT, tt], with the momenta given in units of a^^ . For a 
(110) interface we have instead 



^(k) = -4tcos(fc:r/\/2) cos(fcj^/%/2) - /i. 



(3) 



contains two neighboring atoms which belong to different 
magnetic sublattices A and B. Unit cells chosen below 
for the (100) and (110) interface are shown in Figs. 1 
and 2, respectively. 
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with k^ e [-VStt, VStt], and ky e [-7r/V2, tt/v^], on 
account of the periodic conditions along the crystal sur- 
face. 

The Hamiltonian (1) is quadratic in the Fermi fields 
and can be diagonalized with Bogoliubov transforma- 



tions, c. 



= En("«<T(*)7L + <^V7iaii)%a)- The corre- 
sponding Bogoliubov-de Gennes equations take the form 



E 



V, 



— E„ 



Una{i) 

(4) 

where cr = ± 1 



Here JCf^^^ = -tSi^^-j) + (h^ - n)Sij ± arriiS., 
for up and down spin, Sij and 6(^ij^ are the Kronecker 
delta symbol connecting onsitc and nearest neighbor 
sites, respectively. The off-diagonal block I'ij,^ connects 
the nearest neighbor links P^.o- = — Ay J^y^ with minus 
(plus) signs on the a(6)-links for the c?2,2_j^2-wave pairing 
symmetry, or on-site coupling P^.o- = — A^Jy for conven- 
tional s-wave pairing. We note that a modified Bogoli- 
ubov transformation c^^ = I]„(M^o-(«)7i<^ + Vna{i)%a), 
implemented in Ref. 25, led to modified amplitudes 
Vnaii)' Vna{i) crVnaii)- The Corresponding basic equa- 
tions for these modified amplitudes coincide with those in 
the present paper (in particular, see Eqs.(4), (6)-(9), (13), 
(14)) after redefining the off-diagonal blocks 2?y> — > 



-trl'ij.cr, or, equivalently. A.; 



-crAij, Ai 



Crystal periodicity along the interface makes it conve- 
nient to Fourier transform the BdG equations along the 
y axis and introduce a wave vector component ky. In 
the presence of antiferromagnetic ordering, this should 
be done by taking into account magnetic crystal sym- 
metry along the boundary. In particular, the magnetic 
order parameter nij oscillates rapidly along the (100) in- 
terface and results in a doubling of the period along the 
y-axis. In general, a modified magnetic period along the 
boundary arises for all interface orientations except for 
a (110) interface. We will take into account the mag- 
netic crystal symmetry by introducing a unit cell which 



FIG. 1: (Color online) {100} interface, showing the corre- 
sponding unit cells with two atoms, and basis vectors ei, 62 
of the magnetic lattice. 

On account of the magnetic crystal symmetry, the 
Fourier transformation is taken to be of the form 



A{B) 



A{B) 



27r 



/ MB) 



dky 



Jy 



TXjdy 



{ky) 



Here ky is measured in units of 



a and dy = 2. 



(5) 
V2 

for (100), (110) interfaces, respectively. The transforma- 
tion is identical for atoms of sublattices A and B in the 
same unit cell j. The vector j = {jx,jy) denotes cell 
coordinates, where jx{y) is the x{y) coordinate of the cell 
measured in units of the appropriate basis vectors. For 
definiteness, we identify cell positions with positions of 
the associated site A. 
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cell 



FIG. 2: (Color online) (110) interface, showing the corre- 
sponding unit cells and basis vectors. 

Let vector i denote the location of a nearest neighbor 



4 



site for site A{B) in the same unit cell. Then, the posi- 
tions of all nearest neighbors are described in the case of 
a (100) interface in terms of the basis vectors shown in 
Fig. 1 as (i + ei ± i — ei ± i ± 62, i). Taking this 



into account when performing the Fourier transformation 
(5) in Eqs.(4), we obtain the following one-dimensional 
Bogoliubov-de Gennes equations for the (100) case: 



- ^uf}f\ky) - te±-« {ufiZiky) + uf^fliky) + 2cosfc,<(^)(fc,)) + amf^^^uf}f\ky) - ^f^^^'ik^ 

- e^"" (^S^vfl^M + ^l^^vf^l(k^ + 2cos^,A:^if )^,j^^)(fc,)) = E^^^^ ik^ , (6) 



,vff\k,) + te±-« {vflt^M + -f^tkky) + 2cosfc,.fr)(fc,)) + ^^.^If ^(fc,) - Aff^'u^^iky) 

- e^"" i^^Sir-fltiky) + ^^fn^^M + 2 cos fc.A^^f (fc,)) = Evff^ik^ . (7) 

I 



Here, j denotes only the x coordinate of a cell. The 
factors exp(±ifcy) arise in the nonlocal terms in Eqs.(6), 
(7) since for the (100) interface orientation, in accordance 
with the definitions above, the y coordinate of the site B 
is always less by a than the coordinate of site A in the 
same cell. 

For the (110) orientation, however, the y coordinate of 
the site B is less by a/\/2 than the y coordinate of site A 



in the same cell. Furthermore, if the vector i denotes the 
location of a nearest neighbor site for site A{B) in the 
same unit cell, the positions of all nearest neighbors in the 
case of the (110) interface are described in terms of basis 
vectors shown in Fig. 2 as (i =F Ci, i ± 62, i, i =F Ci ± 62). 
Thus, for the (110) interface orientation we obtain the 
following one-dimensional BdG equations: 



,uff\k,)^2tcoAe ^/2 (uff\k,) + u%t(k,)) + amf^""^ uff\k,) ^ Aff\ff\k,) 



V2 

\A{B).a ±^/2iky , A A(B),b \ ^B{A) 



<r(M+ A- 



A(B),b±^iky 



.A[B).a \ ^B{A) 



EuZ \ky) , (8) 



^.vf^f\k,) + 2t cos ^e^Vl (.g^)(fc,) + v^Uh 



V2 

.A{B).a' ±^/2iky , AA{B),b' 



A(B)A{B} 



iky) 



-Tiky) + 



A^'^hf e^^'ky + A^[^hf) uf^tliky) 



EvTiky) 



(9) 



The singlet superconducting order parameters 



entering Eqs.(6)-(9), are defined as A 



A{B} 



^y^l^^^.A(B),B(A) 

^A(B) 



parameter is m 



.A(B) B(A\ 



The magnetic order 



^MB)^ 



.MB) 



For 



the study of proximity effects, it is convenient to intro- 



duce the magnetization Mi and the pairing amplitude 
Fij, which are related to rrii and A^j by rrii = UiMi and 
Ay ^ -ViF^j, respectively. 

The self-consistency equations in the sublattice repre- 
sentation take the form 



AiB) _ 



E [l<El(fcy)lV(i?«.fc„..) + |<£l(fc,)|V(-£^n,fc.,.) 

n.k-u 



(10) 
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(11) 
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The sum is taken over eigenstates of Eqs.(6)-(9), which 
depend onky, a and possibly an additional set of quan- 
tum numbers n. Eq.(ll) applies to the d-wavc case, 
whereas for the s-wave superconductor with on-site pair- 
ing one should put in Eq.(ll) A"- with amplitudes 

ix o-(^y)' ix si^y) taken for one sublattice a = A,B. 
As usual, f{E) = [1 + exp(£;/T)]~i denotes the Fermi 
distribution function at temperature T . 

Obviously, any bond between nearest neighbors con- 
nects two sites from different sublattices. The notation 
A^ {^fj) means that the order parameter is taken on the 
bond which connects a site of sublattice A (B) within the 
unit cell i with a nearest neighbor site (of sublattice B 
{A)) in the unit cell j. All order parameters are pre- 
sumed to be identical on links (or sites) which can be 
obtained from each other by magnetic translations along 
the interface or interlayer. For this reason, it is sufficient 
to consider i and j in the notation A^^^^ as containing 
only X components, if one indicates in addition the type 
of Hnk (a or b). 



III. ANDREEV EQUATIONS AND 5-MATRICES 

We will base our numerical calculations on the one- 
dimensional Bogoliubov-de Gennes equations formulated 
above, as well as on the corresponding self-consistency 
equations. For an analytical study of superconducting 
interfaces and junctions involving itinerant antiferromag- 
nets, we present in this section the quasiclassical ap- 
proach developed in Ref. 25. As is well known, the 
quasiclassical theory describes physical quantities which 
vary slowly in space compared to the atomic scale, and 
assumes characteristic energies to be much less than 
the Fermi energy Ep. We consider below two types of 
superconductor-antiferromagnet hybrid systems to which 
the quasiclassical approach applies to a different extent. 

The first type of system satisfies the conditions |A| <C 
Ep^ \m\ <^ Ep. The latter inequality guarantees that the 
antiferromagnetic order parameter m, as a rule, varies 
slowly within each separate sublattice A and B. Hence, 
one can use quasiclassical equations both for supercon- 
ducting and antiferromagnetic phases and match them 
at the interface if they are formulated separately for 
each sublattice. Sublattice equations are coupled with 
each other via nonlocal terms which contain, for exam- 
ple, hopping matrix elements or d-wave superconducting 
order parameter fields. In total, this gives us twice the 
usual number of coupled quasiclassical equations. An- 



other possible formulation of the quasiclassical approach 
to the first type of superconductor - itinerant antiferro- 
magnet hybrid systems is not based on the sublattice 
representation. Instead, one can formulate equations 
for quasiparticle trajectories, taking into account that 
the rapidly oscillating antiferromagnetic order parame- 
ter rrij = {—iy'^'^^'"m = exjp(iQj)m couples equations 
for two trajectories, one for a quasiparticle momentum 
kp and the other for {kp + Q). The approach based 
on this {kp,kp 4- Q) representation results also in twice 
the number of quasiclassical equations. There are no fur- 
ther coupled trajectories, since 2Q is assumed to be the 
reciprocal vector of the nonmagnetic crystal due to the 
nesting condition. 

In the second type of hybrid superconductor- 
antiferromagnet systems, only the superconducting order 
parameter satisfies the condition |A| ^ Ep and can be 
safely described with the quasiclassical equations. The 
antiferromagnetic order parameter is taken sufficiently 
large A ^ |m| < £'i? (one could also assume m > Ep 
within the framework of the approach, if this were rel- 
evant). Then the effect of the antiferromagnet on the 
superconductor can be taken into account entirely via 
modified boundary conditions, which complement quasi- 
classical equations at abrupt superconductor - antiferro- 
magnetic interfaces. 



A. Andreev equations in the sublattice 
representation 

We begin with the derivation of Andreev equations in 
the sublattice representation. Assume for this purpose 
that the solution of the BdG equations (6)- (7) and (8)- 
(9) can be represented as the following product of rapidly 
oscillating exponents and a slowly varying Andreev am- 
plitude: 



B 
A 



exp 



kpb 



\ 3a 



\ 



exp{ikp^j:da:j) 



(12) 

Here, kp^^ is the x component of the quasiparticle mo- 
mentum on the Fermi surface, measured in units of a~^. 
The quantity kp^x can be considered a function of ky. 
The Pauli matrix 73 operates in {AS} sublattice space, 
dx = l for (100) interface and = \/2 for the (110) case. 
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Introducing the unit vector b along the crystal 6-axis per- 
mits us to define the Andreev amplitudes in Eq.(12) in a 
unified form which applies to all interface orientations. 

As mentioned above, the parameter fi is considered 
throughout the paper to be small /i <C Ep, since the an- 
tiferromagnetic phase is stable only close to half-filling. 
For this reason, one can additionally include effects 
of small deviations from half-filling in the quasiclassi- 
cal approximation. Taking this into account, we de- 
fine a small parameter of the quasiclassical expansion 
as a = max(|TO|, |A|, |/i|)a/|?Ji?^2:| ^ 1- Here vp^x is the 
a;-component of the Fermi velocity in the normal metal 
state at half filling. The quasiclassical approach works 
only for those ky for which vp^x is not too small and does 
not break the condition a ^ 1. We expand all properties 
associated with the Fermi surface in powers of the small 
parameter {fia/\vp^x\) <SC 1- This concerns, in particular, 
kp,x which enters the exponential factor in Eq.(12). The 
surface adapted Brillouin zone and respective Fermi sur- 
face for small and vanishing /i are shown in Figs. 3 and 
4 for the (100) and (110) interface orientations, respec- 
tively. 




-■k/2 

FIG. 3: (Color online) The Brillouin zone adapted to the (100) 
interface and the Fermi surfaces for small and vanishing jj,. 

For the (100) orientation, the normal-state electron 
band is described as (k) = ^2t(cos fc^-|-cos fcj,) — /i, and 
the Brillouin zone in the case of two atoms in unit cell 
(see Fig. 1) is spanned by fc^; G [— tt, tt], ky £ [— 7r/2, 7r/2]. 
For the half-filled band kp^x = ±(7'' — 1^^!) and Vp^ = 
±2ta sin kp^x- As is seen in Fig. 3, four possible values 
of kp^x, which occur for a given ky ai n ^ 0, merge into 
two values at /i = 0. 



For the (110) interface, we have in the sublattice rep- 
resentation ^*(fe) = =F4t cos(A:j;/\/2) cos(fc,y/\/2) — /i and 

kx.y & [— 7i'/\/2, tt/v^]. If /i > 0, only the Fermi surface 
^~ {k) = exists in the first Brillouin zone, whereas for 
/i < the Fermi surface is determined from ^^(fc) = 0. 
At /i = both Fermi surfaces coincide with the edges 
of the first Brillouin zone. As seen in Fig. 4, two dif- 
ferent values of kp^x, which occur within the first Bril- 
louin zone for a given /cy, touch the edges of the zone 
kp,x — ±n/\/2 in the case of half filling and hence, be- 
come equivalent at ^ = 0. Since the Fermi velocities 
vp^xi and vf^x3 have opposite signs, for the half- filled 
band Vp^{kp) = ±2y/2ta cos (^ky / ^/2) , where two signs 
of vp,x{kp) at the same kp correspond to two degenerate 
parts of the Fermi surface. 
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FIG. 4: (Color online) The Brillouin zone adapted to (110) 
interface and the Fermi surfaces for small and vanishing jj,. 



To obtain the Andreev equations for the (100) orienta- 
tion, we substitute the ansatz (12) into Eqs.(6), (7), sep- 
arately for each of the two values kp,x = ±(7r — ky), and 
disregard terms < a^. For the (110) interface one should 
proceed analogously, using in Eqs.(8), (9) the ansatz (12) 

with kpx ^ TT/V2. We note that - uf^^^l < 

l<-^e+<-^ii-2-2^^l ^"^l-.T^I- Neglect- 
ing terms of the order of , one gets with the required 
accuracy u'^^-^ - = - u'^_-^ = ("j+i ^ "j-i) = 
dxdu'j/dxj. 

Using the outlined procedure, we obtain the following 
Andreev equations in the sublattice representation 



^~B{A) 

,,f,^(S) + I''' + AiB) -.A{B) _ .A{B)' ~A{B) _ l.A{BY , . -BiA) _ ~A{B) , 
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These equations take a unified form, which appfies to 
any interface orientation. For the (100) orientation 

A^^f\kF,x,ky) = 2A^^f\coskF,x - cos ky). The on- 
site d-wave order parameter A^^^\ slowly varying on 
the atomic scale with coordinate j along x axis, is de- 
fined in the coordinate space by the four surrounding 
links. With standard site coordinates i = {ixiiy) one 
can write A^^^ = \ {AaAi+a + \i,ii-a- ^d,ii+b- ^d,ii-b)- 
For the (110) orientation we have A^'^^\kp^xiky) = 

-AA^^f^ sm{kF,x/V2) sin(fc.y/^/2) with kF,x='n-/V2. 

As one can see from the derivation of Eqs.(13), (14), 
they apply to various cases when the sublattice magnetic 
and/or superconducting order parameters vary slowly in 
space, satisfying standard quasiclassical conditions. For 
instance, no particular relation between the sublattice 
magnetizations and is implied yet. For both in- 
terface orientations, the Fermi velocity Vp^{kF) is pos- 
itive for kF,x > within the Brillouin zone in the sub- 
lattice representation. However, the associated solutions 
of Eqs.(13), (14) can describe, in general, both incoming 
and outgoing quasiparticles on either side of the inter- 
face. This is seen from the expression for the density 
of a quasiparticle probability current jp, which can be 
found from the BdG equations in much the same stan- 
dard way known in the case of the Schrodinger equation. 
The probability current density along the x axis, carried 
by the solution with quantum numbers (n, cr, ky), can be 
written in the (100) case as 

'"Fx 

(15) 



Here, the sum is taken over sublattice index v = A, B, as 
well as over the two parts of the Fermi surface a = ±1. 

For the (110) orientation we find 



As usual, the components u and v of the Andreev am- 
plitudes with the same wave vector have opposite contri- 
butions to the probability current'^^. Since the current 
is formed mainly by hopping between nearest neighbor 
sites, it is determined in the sublattice representation as 
a mixed product of A and B components of Andreev am- 
plitudes for any interface orientation. Hence, the sign of 
jp.x depends not only on the crystal wave vector kp, but 
also on the relative signs of the A and B components of 
the Andreev amplitudes. 

One can further transform the equations (13), (14), 
which are formulated in the {j4i?}-sublattice represen- 
tation with two atoms per unit cell into the representa- 
tion with one atom per unit cell. Consider, for example, 
the (110) interface orientation. For one atom per unit 
cell, the Brillouin zone is spanned by kx € [—\/2tt, v^tt], 
ky e [-tt/V^,tt/V2] and ^(fe) is given in Eq.(3), whereas 
in the case of two atoms per unit cell the Brillouin zone 
in the kx direction is: kx G [— 7r/V2, 7r/\/2]. Assum- 
ing \kp_x\ < tt/v^, one can write the following relation 
between the quasiparticle amplitudes in the two repre- 
sentations (compare with Eq.(12)): 



■U2j(fcF + Q) 

U2j+i{kF + Q)e-*('=!'-'^-^.-)/^^. 

r 



(17) 



Here, we have taken into account that if j is the x coor- 
dinate of a two-atom unit cell, then in the representation 
with one atom per unit cell, the site A has even x coor- 
dinate 2j, whereas site B has odd coordinate 2j + 1. 

In the sublattice representation, the wave vector Q = 
(zbv^TT, 0), which is the wave vector of the antiferromag- 
netic pattern that we will study below, is the reciprocal 
crystal vector. Thus, the wave vectors kp and kp + Q are 
equivalent to each other in the approach with two atoms 
per unit cell. In the representation with one atom per 
unit cell, the wave vectors kp and kp + Q are physically 
different. The quantity kp + Q in Eq.(17) is assumed to 



lie in the first Brillouin zone of the representation with 
one atom per unit cell, so that Q = {±^/2tt, 0) should be 
taken there with minus sign for < kp,x < 7r/\/2 and 
with plus sign for —tt/^/2 < kp^x < 0. Thus, it follows 
from Eq. (17) that 

/u/(fcF)\ / U2jikp)+U2jikp+Q) A -."^-^.^ 
\uf{kp)J \u2j+iikp)-U2j+i{kp+Q)J 

(18) 

For the (100) interface orientation, the Brillouin zone 
for the square lattice with one atom per unit cell is 
spanned by kx,y G [—^/2tt, ^/2tt] and ^(fc) is given in 
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Eq.(2), whereas for two atoms per unit cell, shown in 
Fig. 2, the Brillouin zone in the ky direction is: ky G 
[— 7r/-\/2, 7r/-\/2]- Assuming < 7r/-\/2, we find the 

following relation between the quasiparticle amplitudes 
in the two representations: 

[uf{ky)\ ^ 1 .,kj2f^j{ky) + ^jih + Qy)\ .,q^ 

\uf{ky)j 2 \u,{ky)~U,iky+Qy)j ^ ^ ^ 

Here, the wave vector Q = (±7r, ±7r) is the reciprocal 
crystal vector in the sublattice representation. The quan- 



tity ky + Qy in Eq.(19) is assumed to lie in the first 
Brillouin zone of the representation with one atom per 
unit cell, so that Qy = ±7r should be taken with mi- 
nus sign for < kp.y < 7r/\/2 and with plus sign for 
-n/V2 < kp.y < 0. 

Substituting Eq.(18) or (19) into Eqs. (13), (14), 
we obtain the Andreev equations in the {kpjkp + Q)- 
representation. Since quasiparticle crystal momenta kp 
and kp+Q both lie on the Fermi surfaces on either side of 
the interface, the Andreev equation in the [kp^kp + Q)- 
representation takes the comparatively simple form: 



(^-iVp,x{kp)p^Tz-^ - ^JLT^ 



A R 

rUj — m 



^h{kp)—pM,,^E,i>,, 



(20) 



Here, {ujcr{kp),Uja{kp -f Q),Vja{kp),Vja{kp + 

Q)). Ta and Pa denote the Pauli matrices in particle- 
hole and {kp, kp + Q} spaces. Eq.(20) applies generally 
for any particular relation between m"^ and 'm^ ■ For 
instance, it applies also to the study of weak ferromag- 
nets (m <g; e/). In the case of antiferromagnetic ordering 
satisfying the condition = — w^, within the quasi- 
classical accuracy one can disregard in Eq.(20) the term 
containing mf + mf . 

The only term which couples Andreev amplitudes with 
momenta kp and kp + Q in Eq.(20) contains the differ- 
ence between sublattice magnetizations mf — mf . This 
is natural, since a finite difference 7nf — mf results in a 
doubling of the period in the system which we study. Eqs. 
(13), (14) and (20) apply also in the absence of period 
doubling, being equivalent in this case to the standard 
Andreev equations. If period doubling takes place only 
at the boundaries, the sublattice or {kp,kp + Q) rep- 
resentations of the quasiclassical equations can be con- 
venient for applying appropriate boundary conditions to 
the solutions. 



interface can be written in the following form 



B. Boundary conditions and S-matrix in 
{kpjkp + Q) representation 

The assumption of slowly varying order parameters 
does not apply in the vicinity of abrupt boundaries. This 
invalidates quasiclassical equations close to the bound- 
aries and makes it necessary to complement them with 
appropriate boundary conditions. The conditions have 
been obtained for Andreev amplitudes in Ref. 33 and re- 
derived later by various methods (see, e.g. Ref. 34). The 
boundary conditions for Andreev amplitudes at a plane 



*1 



5*11 5i2 

'5'21 S22 



(21) 



Here, ^' denotes a collection of Andreev am- 
plitudes. For example, in the {kp,kp + Q) 
representation it contains eight amplitudes 
= {uja{kp),Uj„{kp + Q),Vjs{kp),Vjs{kp + 
Q),Uja{kp),Uja{kp + Q),Vjcr{kp),Vja-{kp + Q)) . The 
superscripts l{r) indicate that the amplitudes are taken 
on the left (right) side of the boundary. The subscripts 
± in the amplitudes denote the sign of the Fermi velocity 
components vp,:c{kp) or vp,x{kp -\- Q) for electrons. 
Thus, the solutions entering the left and right hand sides 
of Eq.(21) are connected at the interface by the normal- 
state scattering S matrix: S = \\Sij\\ = 1,2). This 
matrix Su contains the reflection amplitudes of normal- 
state quasiparticles from the interface in i-th half-space, 
whereas Sij with i j incorporates the transmission 
amplitudes of normal-state quasiparticles from side j. 
In the (kp,kp -\- Q) representation, each component Sij 
= 1,2) in the matrix 5 is an 8 x 8 matrix in the 
eight-dimensional product space of particle-hole, spin 
and {kp,kp + Q) variables. We introduce the Pauli 
matrices Pa-, and ctq, which operate in {kp,kp + Q} 
space, particle-hole space and spin space respectively. 

The normal-state S matrix is diagonal in particle-hole 
space: 



5=1 f^]^s^-±ll + S- 
021 022 



(22) 



If the AF order parameter does not change its direc- 
tion, one can take a quantization axis along m. In this 
case up and down spin states are eigenstates of the BdG 
and Andreev equations, which are formulated separately 
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for each electron spin orientation as given in Eqs.(4), (6)- 
(9), (13), (14). Then the associated S matrix turns out 
to be diagonal also in spin space, and the boundary con- 
ditions reduce to the following equalities: 

(23) 

~a.r f ca/3 -3.1 , nap ~J3,r \ 

(24) 

K 13 ^ ' 

Here, the superscripts a, /3 take the two values kp.yi 

kF,y + Qy ■ 

As an example, we apply in the following the quasiclas- 
sical approach to analytical calculations of the subgap 
spectrum of quasiparticle interface states near a (110) 
SC/AF interface in the absence of potential barriers. A 
detailed self-consistent study of quasiparticle states at 
interfaces with antiferromagnets will be presented be- 
low in Sec. IV and V. Let an s-wave superconductor 
and an antiferromagnet, separated with the (110) inter- 
face, occupy separately the right {j > 0) and the left 
(j ^ 0) half-spaces of the square lattice. Assuming 
\E\ < |A| ^ \m\, \vF,x\/a, we use the Andreev equations 
only in the superconducting region and apply the bound- 
ary conditions at the superconductor-antiferromagnct in- 
terface. Wave functions for quasiparticlcs with energies 
below the superconducting gap decay exponentially with 
increasing distance from the interface. Solving the An- 
dreev equations, one can easily find the standard two- 
component solutions ip'^^j.{kp) = {u]j. j.{kp),Vg j.{kp)^ 
taken in the right half-space of the boundary: 

«,(M^C,(^±'^J^)^ ,25, 

Reflection amplitudes for electrons r'^ and holes rj, 
which enter the quasiclassical boundary conditions, are 
taken for the normal-metal state of the superconducting 
region at the Fermi surface. An AF/N boundary is im- 
penetrable for normal-metal quasiparticlcs with energies 
below the antiferromagnetic gap, even in the absence of 
any interface potentials (i.e. for a transparent interface). 
Hence, the corresponding transmission amplitudes van- 
ish and the complex reflection amplitudes for electrons 
and holes have unit modulus: |r'^(fcF)| = \r'^{kp)\ = 1. 
Further, if ip^{kF) represents an outgoing solution in the 
case of a (110) interface, then 'tj}^{kp + Q) is the incoming 
solution, and vice versa. Here, Q = (±V27r, 0) only has 
nonzero x component. For this reason, it is convenient to 
formulate the boundary conditions for the (110) interface 



indicating explicitly only the ky component: 

(26) 

A relation between r^{kp) and r'^{kp) follows from the 
fact that the quantity —[u-cr{kp)]* represents the wave 
function v^i^kp) for a hole with energy —E from one spin 
sub-band if u^{kp) is the wave function for an electron 
with energy E from the other spin sub-band. Under the 
condition A ^ t, one can consider the reflection am- 
plitudes at subgap energies ±£' as taken at the Fermi 
surface. Accounting additionally for the fact that in- 
coming and outgoing waves for holes are interchanged as 
compared with the case of electrons, we find for normal- 
metal electrons and holes at the (110) AF/N interface: 

= (r| ) . This condition simplifies because of the 
equalities |r'^(fci?)| = \r'''{kp)\ = 1: = ^-a- Eventu- 
ally, applying the boundary conditions (26) to the solu- 
tion (25), we get 

C+ (e+ ^v/A^~P) = rUkp)C^ (e- ^v/A^~P) 

(27) 

This equation determines the energies of the bound states 
at a (110) AF/sSC interface. In order to present the en- 
ergy in a convenient form, we divide the reflection am- 
plitude into two parts r^{kp) = r^p + rg^^, which are 
symmetric and antisymmetric with respect to spin in- 
version, respectively. The first part r^p is actually spin- 
independent and can be considered as the contribution 
to the reflection amplitude from the specular reflection 
channel. The contribution to the reflection amplitude 
from the spin-dependent Q reflection differs in symme- 
try from the specular reflection part. This is the part 
antisymmetric in tr, rg^o-, which is an imaginary quan- 
tity. The phase of rq^ differs by tt for spin up and 
down quasiparticlcs^'^. It follows from the definition 
given and the condition |r^(fci?)| = 1 that reflection coef- 
ficients Rsp = l^'spP and Rq = I^'q.ctP satisfy the relation 
Rsp + Rq = 1. With these quantities, we obtain from 
Eq.(27) the following bound state energies at the (110) 
AF/sSC interface: 

E = ±A,yJ R,p{kp) . (28) 

If a d-wave superconductor occupies the right half 
space instead of the s-wave one, the expression for 
ip'^ j.{kp) in Eq.(25) is modified with the substitution 
As ±Ad{kp). As a result, we obtain the spectrum of 
the bound states at the (110) AF/dSC interface: 

E = ±Ad{kp)^RQ{kp) . (29) 

The specific expressions for the normal-state reflection 
amplitudes can be found for the (110) AF/N interface 
along the standard way by solving the Schrodinger equa- 
tions for electrons in the left and the right half-spaces 
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and constructing incoming and outgoing solutions in the 
normal-metal region, as well as exponentially decaying 
solution in the antiferromagnetic region for electrons with 
energies below the antiferromagnetic gap. The reflection 
amplitude is fixed after substituting the bulk solutions to 
the equations, taken at lines where the nearest neighbor 
hopping mixes the normal-metal and the antiferromag- 
netic regions: 



between AF and either s-wave or d-wave superconduc- 
tors. The coupling constants are site dependent with: 



rspikp) = 



-1/2 



y2vF,xiky) ^ 

( V^VF,x{ky) 

\ ma 



-1/2 



(30) 



where VF,x{ky) = 2\/2ta cofi{ky / \/2) . 

Substituting Eq.(30) into Eqs. (28), (29), we come 
to the results obtained in Ref. 25. We note now that 
the (110) interface represents a special situation, when 
the normal state Fermi velocity for the half-filled lat- 
tice has only a vp^^ component, which is perpendicular 
to the surface. Hence, outgoing and incoming normal- 
metal quasiparticles move only along or opposite to the 
surface normal. The quasiparticle trajectory of this kind 
is intrinsic to the specular reflection channel. However, 
Q reflection takes place in this particular case along the 
same quasiparticle trajectory. For this reason, specular 
and Q reflections make coherent contributions to the to- 
tal reflection amplitude rj^(fei?) = Tgp+rQ fj. For different 
interface orientations, for example, for the (100) inter- 
face, specular and Q reflection takes place along different 
trajectories and do not interfere with each other. 

In the absence of interface potentials specular reflec- 
tion of quasiparticles arises entirely due to a mismatch 
of Fermi velocities in the AF and the sSC. Since nor- 
mal metal states are presumably identical in the left and 
right halfspaces under the conditions A <C t the mis- 
match in the model is controlled by the parameter m/t. 
As is seen from Eq.(30), Tsp —>■ 1 when the antiferro- 
magnetic gap is large, \m\ At cos ky. At the same 
time, Q reflection becomes dominant in the opposite limit 
|m| ^ At cos ky, taking place at energies below the an- 
tiferromagnetic gap \E\ < \m\. The quasiparticle bound 
state energy (28) at the AF/sSC interface is almost zero 
in the limit \E\ ^ As, whereas the bound state (29) at 
the AF/dSC interface hes very close to the edge of the 
superconducting gap. 



u 

V 



for 
for 



i < 0, 
i > 0, 



(31) 
(32) 



and zero elsewhere. In the following we will typically be 
studying finite systems of length N = 100 — 200 along 
the X axis. The ky sum is performed explicitly by using 
400 points in the Brillouin zone. The interface is always 
positioned at the bond in the middle of the system and 
the potential hi is only non-zero on the two sites immedi- 
ately adjacent to the interface. We apply open boundary 
conditions equivalent to an impenetrable wall at each end 
of the system. 



s-wave superconductor - antiferromagnet (100) 
interface 



We begin with the AF/sSC (100) situation since Q re- 
flection is expected to lead to low-energy bound states for 
this particular interface orientation^^. As is well known, 
nonmagnetic interfaces do not break s-wave Cooper pairs 
in their vicinity. In contrast, ferromagnetic as well as an- 
tiferromagnetic interfaces are pair breaking, in general, 
even for s-wave superconductors since they break time- 
reversal symmetry. This is analogous to the difference 
between efi^ects of nonmagnetic and magnetic impurities 
in s-wave superconductors. In Fig. 5 we show the self- 
consistent suppression of the magnetization Mi and the 
pairing amplitude Fa near the (100) interface for vari- 
ous values of superconducting and antiferromagnetic cou- 
pling constants. As seen, the healing length on each side 
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FIG. 5: (Color online) Self-consistent spatial dependence of 
the absolute values of the magnetization Mi and the pairing 
amplitude Fa near a (100) interface. Parameters: fi = and, 
from top to bottom, {U = 2.0t,V = 2.0t), {U = 1.6t,V = 
1.5t), {U = 1.2t,V = 1.2t), and (C7 = 0.93t,V = 0.9t). 



IV. QUASIPARTICLE STATES AT SC/AF 
INTERFACES 

In this section, we study quasiparticle spectra and the 
corresponding local density of states in the vicinity of 
AF/SC interfaces based on our self-consistent approach 
outlined above. We consider (100) and (110) interfaces 



of the interface decreases with increasing amplitude of 
the order parameter in agreement with the behavior of 
the respective magnetic and superconducting coherence 
lengths hvF,x/\'m\, hvF,x/^s- We find a correlation be- 
tween the strength of the suppression of order parameters 
and the energy of the Andreev bound state arising near 
the AF/sSC interfaces. The lower the (positive) energy 
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of the subgap state, the stronger the suppression of both 
order parameters at the interface, at least in simple cases 
which have been studied. 



where 




FIG. 6: Eigenvalues for the (100) AF/sSC interface as a func- 
tion of ky for 11 = 0.0, U = 2.7i, Vs = 2.0t and h = 0.0 (a) 
and h — 2.0t (b). Here, one sees explicitly the presence and 
dispersion of the bound state band inside the gap. 



In Fig. 6 we plot the quasiparticle spectrum as ob- 
tained from the eigenvalues of the BdG equations (6)- (7). 
Naturally, bound states present at the interface show up 
inside the main gap of the spectrum as a distinct band, 
which disperses with the momentum component ky along 
the interface. The two gap edges seen in Fig. 6 are asso- 
ciated with the superconducting (lesser) and the antifer- 
romagnetic (larger) gaps. 

We have calculated the bound state spectrum also an- 
alytically, assuming <^ m, t. Solving the Andreev 
equations for the superconducting region and applying 
appropriate boundary conditions, we obtain the two dis- 
persive energies of quasiparticle Andreev bound states for 
the (100) AF/sSC interface. The spectrum is symmet- 
ric with respect to the zero level and can be described 
with a reflection coefficient Rsp{ky) for quasiparticles in 
the specular reflection channel: E{ky) — ±Asy/ Rsp{ky). 
This expression actually coincides with Eq.(28) derived 
in Sec. Ill for the (110) interface although the explicit 
expression for the reflection coefficient differ for the two 
orientations. A calculation of the reflection coefficient 
Rsp{ky) for the (100) AF/sSC interface in the absence of 
potential barriers leads to the following explicit expres- 
sion for energies of the Andreev bound states: 



E{ky) = ±As 

^ A{ky) 



AHky)+4[f^ 



1/2 



A{ky) + 2 sin^ ky + ^A^ikv) + '^{^y J 



(33) 



Aiky) 



\2t 



V 



(34) 



The dispersion shown in Fig. 6a can be verified to agree 
well with the expression in Eq.(33) within the accuracy 
±(A/t)2. Eq.(33) is very similar to Eq.(7) of Ref. 35 
and can be obtained from there simply by substituting 
the magnetic m and s-wave order parameters for the 
charge density wave Ws and d-wave A^ order parameters 
respectively. As follows from Eq.(33), the quasiparticle 
subgap state becomes a dispersionless zero-energy state 
if one additionally assumes m <C i and disregards terms 
less or the order of {m/t)'^. This limiting case corresponds 
to the zero-energy solution found from the quasiclassical 
Andreev equations applied to both superconducting and 
antiferromagnetic regions under the conditions to, A^ <^ 

The differential tunneling conductance measured, for 
instance, by scanning tunneling microscopy (STM) ex- 
periments is proportional to the local density of states'^^. 
Therefore, it is important to calculate the LDOS, given 

by 



Im x - 

TT ^ 

nkyCr 



rikyW 



,iky)\' 



En 



■ En 



-iV 
(35) 

where a = A^ B indicates the magnetic sublattice and F 
is an artificial broadening which in the following is set 
to F 0.02<. For the (100) interface, we find Nf'{uj) = 
Nf{uj). In plots of the resulting LDOS we expect any 
bound states to result in additional peaks inside the gap 
of the bulk AF and SC. These peaks should be localized 
near the interface. This is indeed the case for the 100 
AF/sSC interface, as can be seen from Fig. 7. Here, the 
two center LDOS scans in both (a) and (b) are at the 
interface while the top (bottom) five scans are moving 
into the SC (AF). 

Additional potentials h ^ present near the interface 
can strongly enhance specular reflection at the expense 
of Q-reflcction. There are several consequences of the 
interface potential for Andreev bound states present in 
the system. First, potentials tend to suppress the bound 
states resulting from Q-reflection and move their posi- 
tions towards the gap edge. As expected, in the limit 
h ^ t wc always find that the Q-reflection bound states 
have been pushed into the continuum. Second, in the 
regime where h is of the order of t, we find that the 
main effect of the specular reflection channel is to cause 
a stronger dispersion of the bound state energy. One can 
identify additional extrema in the wave vector depen- 
dence of the bound state energy E{ky). A typical exam- 
ple is seen in Fig. 6b where h = 2.0t. The new stationary 
points in the dispersion lead to additional LDOS peaks 
near the interface as seen in Fig. 7b. In the LDOS we 
also see that the particle-hole symmetry is broken when 
h 0, whereas the quasiparticle spectrum is still sym- 
metric with respect to the Fermi level. A similar asym- 
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FIG. 7: (Color online) LDOS corresponding to the same pa- 
rameters as in Fig. 6. The interface bound states result in 
sub-gap peaks in the LDOS near the interface region. The 
two center LDOS scans in both (a) and (b) are at the inter- 
face while the top (bottom) five scans are into the SC (AF). 
The lines are off-set for clarity. 



mctry between positive and negative bias in the LDOS 
will be present when starting from a particle-hole asym- 
metric band, i.e. when /x 7^ 0. For the sake of clarity, 
the results presented below are for the particle-hole sym- 
metric nested band where fj, = and any asymmetry will 
only result from a non-zero interface potential h. 

It is interesting to investigate the importance of the 
suppression of the order parameters near the interface re- 
sulting from the self-consistency. Fig. 8 shows the bands 
and the corresponding LDOS for a non-self-consistcnt 
calculation with step-function fields: mf ~ ^"m-f = 
77ioO(— ?) and = AoO(i) with mo — 0.7t and Ag = 
OAt. Clearly, the results are very similar to those shown 
for = in Fig. 6a and Fig. 7a, respectively. This 
result applies also to the other interfaces studied below: 
in general the self-consistency has only minor effects on 
the bound states resulting from Q-reflection. It can, on 
the other hand, depending on specific parameters, induce 
new bound states close to the gap edge. 

As a final point in this section, we verify the bound na- 
ture of the subgap states by showing explicitly the spatial 
dependence of the eigenstates corresponding to the sub- 
gap band in e.g. Fig. 6a. This is done in Fig. 9, where 
we plot \uf{ky)\'^ = Wf{ky)\'^ as a function of the x- 
componcnt of the unit cell near the interface for ky = 0.0 
and ky = 0.257r. Clearly, the wavefunctions are seen to 
be bound to the interface region. 
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FIG. 8: (Color online) Non-self-consistent spectrum and 
LDOS corresponding (roughly) to the same parameters as in 
Fig. 6a but with step-function spatial dependence of Mi and 
Fa. These results are basically identical to those shown in 
Fig. 6a verifying that it is the Q-reflection, not the order 
parameter suppression, that generates the subgap states. 
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FIG. 9: Amplitude of |u^(fcj;)p as a function of the x- 
component of the unit cell i corresponding to the bound state 
in Fig. 6a when ky — 0.0 (a) and ky — 0.257r (b), respectively. 



B. s-wave superconductor - antiferromagnet (110) 
interface 

Unlike the (100) interface studied in the preceding sec- 
tion, in the (110) case all sites at the interface belong to 
the same sublattice, so that the interface layer (chain) 
itself is ferromagnetically ordered. Up and down magne- 
tizations alternate only along the interface normal, i.e. 
the x-direction as seen from Fig. 2. For this reason 
some characteristic properties of AF/sSC (110) interfaces 
could naively be expected to be reminiscent of the prop- 
erties of superconductor-ferromagnetic boundaries. For 
example. Cooper pair wave functions are known to de- 
cay into the ferromagnet adjacent to the superconductor, 
manifesting at the same time spatial oscillations^'^^. The 
oscillations are known to be induced by the difference be- 
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tween the momenta of spin up and down quasiparticles 
with the same energy^^. These oscihations in the SC/F 
proximity effect are related to those in the Fulde-FerreU- 
Larkin-Ovchinnikov (FFLO) superconducting phase'^®. 

In Fig. 10 we show the obtained self-consistent spa- 
tial profiles of the magnetization Mi and the pairing am- 
plitude Fa near the (110) AF/sSC interface region for 
various coupling strengths. Both the magnetization and 
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FIG. 10: (Color online) Spatial dependence of the self- 
consistent results for the absolute value of the magnetization 
Mi and pairing amplitude Fu for the AF/sSC 110 interface 
orientation. Parameters; fi = and (top) U = 2.0t, V = 2.0t, 
(middle) U = 1.6t, V = 1.5t, and (bottom) U = 0.93t, 
V = 0.9t. 

the pairing amplitude display an oscillatory decaying be- 
havior near the interface. The oscillations have even-odd 
character, i.e. they are not present within each separate 
magnetic sublattice. Thus, the oscillations in Fig. 10 are 
not equivalent to FFLO oscillations, but simply induced 
by the AF staggered ordering. The characteristic scales 
for suppression of the order parameters are seen to be the 
corresponding coherence lengths of the antiferromagnet 
and the superconductor. We find that an additional po- 
tential at the interface decreases both the decay length 
and the amplitude of the oscillations on both sides of the 
interface. 

The dispersion of the subgap energies at (110) AF/sSC 
interfaces are described analytically by Eqs.(28), (30), 
within the non-self-consistent quasiclassical approach. In 
Fig. 11(a) we plot the numerically determined self- 
consistent (spin up) eigenbands for the (110) AF/sSC 
interface. The dispersion of the main bound state can be 
checked to yield excellent agreement with the result in 
Eq. (28). Positive and negative energies correspond there 
to spin up and down quasiparticles respectively. Thus, at 
low temperatures only spin down quasiparticlc Andrcev 
states will be occupied. Strongly spin-discriminated An- 
dreev states together with the alternating magnetization 
of chains, which are parallel to the (110) interface, are 
at the origin of the even-odd oscillations of the order pa- 
rameters shown in Fig. 10. It is worth noting that spin 
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FIG. 11: Spin up eigenvalues resulting from the Bogoliubov- 
de Gennes equations in the case of a 110 AF/sSC interface. 
The parameters are U = 2.7t, V = 2.0t, and h = 0.0 for figure 
(a) and h = l.Ot for (b). The bands shown here correspond to 
{ui,vi). The bands associated with (m|,D|) are identical to 
the bands shown here upon mirror reflection around E — 0. 

polarized Andreev states are compatible with singlet spin 
structure of Cooper pairs^. Whereas the electron and 
the Andreev reflected hole belong to different spin sub- 
bands and have opposite Zeeman energies, they possess 
identical spin polarization. Since they have also almost 
opposite velocities, they do not carry together any spin 
current except for special cases'*". 

At the edge of the Brillouin zone we find new 
high-energy subgap states. By comparing to non-self- 
consistent calculations with step-like spatial order pa- 
rameter dependence, we have found that these states are 
related to the self-consistent order parameter suppression 
near the interface. When including a potential at the 
interface, the structure of the dispersion of the subgap 
states becomes more complicated. In Fig. lib we show 
an example where h = t. When the potential becomes 
too large, /i ^ m. A, the specular refiection channel com- 
pletely dominates the Q-reflection and all bound states 
arc pushed out of the gap into the continuum. 

The spin-summed LDOS corresponding to the results 
shown in Fig. 11a and Fig. lib is shown in Fig. 12. 
Since in the (110) case sites A and B from the same cell 
are at different distances from the interface, Nf-{uj) and 
A^j^ (a; ) differ from each other. Each of Figs. 11a, lib dis- 
play Nf{Lu) and Nf{uj) jointly, i.e. shows the LDOS for 
all sites along the x-axis. As expected the bound states 
again show up as peaks in the LDOS close to the inter- 
face. The states at the Brillouin zone edge can give rise 
to higher energy peaks. The main influence of the po- 
tential scattering is to break the particle-hole symmetry 
present in the LDOS in Fig. 12a. In particular, as seen 
from Fig. 12b, the presence of a potential h can lead to 
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distinct even/odd amplitude modulations of the bound 
state LDOS peaks into the superconductor. 
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FIG. 12: (Color online) LDOS near the (110) AF/sSC inter- 
face corresponding to the parameters used in Fig. 11. The 
two center LDOS scans in both (a) and (b) are at the interface 
while the top (bottom) five scans are the LDOS upon moving 
into the SC (AF). The lines are off-set for clarity. 



C. d-wave superconductor - antiferromagnet (100) 
interface 

The electronic structure of interfaces formed by d-wave 
superconductors and antiferromagnets is an important 
problem relevant to e.g. the high temperature supercon- 
ductors where AF and dSC order dominate the phase 
diagram. A qualitative difference between the properties 
of AF/sSC and AF/dSC interfaces arises from the fact 
that the quasiclassical c?-wave order parameter A^(fc/) 
changes sign in a Q reflection event for any interface 
orientation: Ad{kp + Q) = —Ad{kp). This change of 
sign can strongly weaken the effect of phase difference 
of reflection coefficients for spin up and down quasipar- 
ticles, which is close to tt in the limit A, to <C For 
this reason, low-energy interface bound states Esikp) ^ 
min{TO, Abuik(feF)} existing at AF/sSC interfaces under 
the conditions Ag, m ^ t, do not exist at AF/dSC in- 
terfaces with arbitrary orientation^^. In the present sec- 
tion we demonstrate results of self-consistent numerical 
calculations for the spatial profiles of superconducting 
and antiferromagnetic order parameters, the quasiparti- 
cle spectrum, and the associated LDOS in the vicinity of 
(100) AF/dSC interfaces. 

Fig. 13 shows a typical result for the suppres- 
sion of Mi and Fi near the interface. Here, we 
plot the d-wave pairing amplitude defined on-site in 
the usual way from the four surrounding links Fi ~ 
\ {Fi^i+a + Fis-a - F,^i+b - F,^b)- Comparing Fig. 13 
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FIG. 13: (Color online) Spatial dependence of the self- 
consistent results for the absolute value of the magnetization 
Mi and pairing amplitude Fi for the AF/dSC 100 orientation. 
Parameters: = 0, and (top to bottom) U = 2.7t, V — 2.0t; 
U = 1.7t, V = l.Ot; and U = l.lt, V = 0.5t. 



and Fig. 5, one can notice weaker pair breaking effects 
on the scale of the respective coherence length for the 
AF/dSC (100) interface compared to the AF/sSC inter- 
face. This is associated with lower bound state ener- 
gies at AF/sSC interfaces. We find similar results for 
the (110) AF/dSC interface studied in the next section. 
Contrary to the (110) AF/sSC case, the order parame- 
ters exhibit a rather smooth suppression near a (100) or 
(110) AF/dSC interface, i.e. the order parameter oscil- 
lations are absent. This correlates with the absence of a 
spin discrimination in the quasiparticle subgap spectrum 
generated by the (110) AF/dSC boundary. 

The eigenbands for two different parameter sets are 
shown in Fig. 14. The left figure shows the result without 
any potential at the interface. Here, a subgap band exists 
close to the edge of the continuum. We find that the 
larger the ratio m/A, the closer the subgap band is to 
the gap edge. This agrees with Ref. 25, where it was 
demonstrated that there are no interface bound states 
on AF/dSC (100) interfaces in the limit to > A. Only 
upon decreasing the ratio m/A does a bound state gets 
peeled off the continuum, as seen in Fig. 14a. The band 
is degenerate and will be split by a non-zero potential 
/i 7^ as seen in Fig. 14b. This is accompanied by the 
appearance of additional extrema in the dispersion of a 
lower band. Although with increasing h a lower band 
first becomes closer to the Fermi level than for ft, = 0, 
its position depends nonmonotonically on h and never 
approaches zero energy. When h ^ t both bound state 
bands are pushed into the continuum. The results of non- 
self-consistent calculations are very similar to the results 
shown in Fig. 14: it is the Q-reflection channel that 
induces the bound states at the AF/dSC interface, not 
the order parameter suppression. 

The LDOS corresponding the same parameters used in 
Fig. 14 is shown in Fig. 15. Again the new LDOS peaks 
arising from the bound states are sensitive to the presence 
of the interface potential h. When h = the bands in 
Fig. 14a close to the continuum edge gives rise to LDOS 
peaks near the coherence peaks. Experimentally, it may 
be a challenge to distinguish these bound state peaks 
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(a) (b) 




FIG. 14: Eigenbands for the AF/dSC 100 interface as a func- 
tion of ky. Parameters: U = 2.7t, V = 2.0t, jj. — and 
h = 0.0 (a) and h = l.Ot (right). 
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FIG. 15: (Color onhne) LDOS corresponding to the plots 
shown in Fig. 14. The two center LDOS scans in both (a) 
and (b) are at the interface while the top (bottom) five scans 
shows the LDOS when moving into the SC (AF). 

from the coherence peaks of the bulk dSC. However, for 
the case when h — l.Ot, we see that the lower band from 
Fig. 14b results in sharp LDOS peaks near the interface 
in the intermediate region of subgap energies. 



D. d-wave superconductor - antiferromagnet (110) 
interface 

As is well-known, a zero-energy Andreev bound state 
exists at the (110) insulator I/dSC interface generated by 



the sign-reversal of the gap function as seen by a quasi- 
particle specularly scattered off the surface^^"'^^. This 
state has been observed in the differential tunneling con- 
ductance as a conductance peak at zero bias^^"^^. In 
this section we study the AF/dSC interface with (110) 
orientation. 

In Fig. 16 we plot the (spin up) bands and the corre- 
sponding LDOS when U = 2.7t, V = 2.0t and h = 0.0. 
Since the open boundary conditions at the edges of our 
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FIG. 16: (Color online) Bands and LDOS for the 110 AF/dSC 
interface with U = 2.7t, V = 2.0t and h = 0.0. The low- 
energy bands are dominated by the gap node at ky = 0, the 
Q-reflection bounds states and a zero energy state. The lat- 
ter state is located at the dSC/I boundary at the right-most 
the end of our system where the open boundary conditions 
operate as a hard wall. 

finite system are equivalent to a hard wall, a zero en- 
ergy state (ZES) is present at the superconducting end 
of our system. This state, clearly seen in Fig. 16a, is 
not associated with the AF/dSC interface in which we 
are interested. It mixes with the bound states result- 
ing from the Q-reflection at the interface, however, and 
causes small deviations from the expected mirror sym- 
metry (through = 0) of the bound state bands. In the 
limit of an infinite system, this spurious effect disappears 
and the continuum gap node closes at ky — 0. 

The dispersion of the subgap states in the ab- 
sence of potential barriers is described analytically in 
Eqs.(29), (30) assuming spatially constant order param- 
eters Ad{kp) <C TO, t. These states should move to 
lower energy as to increases. For our band vf,x = 
2y/2ta coa{ky / ^/2) and again it can be easily verified that 
the functional form of the self-consistent subgap band in 
Fig. 16a agrees well with Eqs. (29), (30). 

The presence of a ZES at (110) I/dSC interfaces, makes 
it interesting to plot the LDOS as a function of increased 
interface potential h. Fig. 17 shows the evolution of the 
LDOS when approaching the interface in two cases where 
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FIG. 17: (Color online) LDOS for the 110 AF/dSC interface 
with fi = 0,U = 2.7t, V = 2.0t and h = l.Ot (a) and h = 2.0t 
(b). Again, the two center LDOS scans in both (a) and (b) 
are at the interface while the top (bottom) five scans shows 
the LDOS when moving into the SC (AF). Here, one clearly 
sees the continuous evolution of the well-known zero energy 
state with increasing interface potential h. 



h = t (a) and h = 2t (b). The details of the dispersion of 
the subgap states are sensitive to the strength of the po- 
tential barrier at the interface, and the resulting LDOS 
will strongly depend on h. In the limit where h ^ t the 
low-energy LDOS near the interface will be dominated by 
the ZES. However, in the regime where /i ~ i, the ZES 
coexists with the Q-reflection bomrd states as is evident 
from Fig. 17. This finding is relevant for the discussion 
of possible surface induced magnetization near I/dSC in- 
terfaces in cupratc superconductors^^'^^'^^. In the case 
of finite Mi near a (110) dSC surface, we would expect 
small side-band peaks originating from the Q-reflection 
as seen in Fig. 17. To the best of our knowledge no such 
features have yet been identified in experiments. 



E. Transfer matrix method 

The existence and dispersion of bound states at inter- 
faces between SC and AF can also be conveniently formu- 
lated within a transfer matrix method designed to locate 
the bound states from their defining property: spatially 
decaying wavcfunctions. Below, we use the same transfer 
matrix formalism presented in Ref. 24. In this method, 
one introduces a (g, e)-dependent matrix T {i + l,i) de- 
fined by 



(36) 



which transfers the spinor ^ from site i to site i ~\- 1. 
For a model with nearest neighbor coupling ^' takes the 



explicit form 5* (i) = (■(/; (i) ,ip{i — 1)) where 

■0 (i) = {Ukya {i) , Vk^a (i) , Uky+yra («) , {i)) . (37) 

The associated 8x8 transfer matrix has the general form 



T{i + l,t) 



A B 
1 



(38) 



where A (B) denotes the 4x4 coefficient-matrix con- 
necting ip {i + 1) and "0 (i) {ip (i — 1)) determined from 
the Bogoliubov-de Gennes equations (4). 

For a typical interface there will be three distinct trans- 
fer matrices; one in the bulk magnetic region Tm, one in 
the bulk SC region Tsc a-nd one associated with trans- 
fer through the interface T/. By diagonalizing Tm and 
Tsc one determines whether eigenstates decay, grow or 
propagate from the interface depending on whether the 
eigenvalues are less than, larger or equal to one, respec- 
tively. Here, decaying and growing refer to propagation 
along the x-axis. 

Let PEm denote the matrix obtained after propagat- 
ing the eigenvectors of the bulk magnetic transfer matrix 
through the interface. Then wc can define a matrix A 
given by 



PE 



M 



Esc- A 



(39) 



where Esc is the matrix containing the eigenvectors of 
the bulk superconducting region as column vectors and 
dot indicates matrix multiplication. Let Sg"^ and 5'^'-^ 
denote the subspace of growing eigenstates of PEm and 
Esc respectively. Consider the following linear combina- 
tion of the growing states of PEm 

Y,P,\PEMi> = P^MEsc]> (40) 



J2 E ^^-^^^ i^^cj > 



^From equation (40), it is evident that to have a bound 
state at the interface, the vector (3 must belong to the 
null space of the reduced matrix Ar, defined to be the 
Sg'^ X Sg^ upper left part of the original matrix A. This 
follows since the matrices PEm and Esc were chosen 
to have the eigenstates with the largest eigenvalues as 
column vectors to the left. Thus, when the two subspaces 
Sg'~' and Sg^ have the same dimension, a bound state 
at the interface is characterized by the vanishing of the 
determinant of Ar 



Bound state criterion: 



det Ar = 0. 



(41) 



This is the criterion used previously by Andersen and 
Hedegard to study the splitting of the ZES near dSC/AF 
interfaces^'*. 

In Fig. 18, we show the results of Eq. (41) for the 100 
(a) and 110 (b) AF/sSC interface, respectively. Here, we 
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FIG. 18: Comparison of the dispersion obtained for the 100 
(a) and 110 (b) AF/sSC interface form the transfer matrix 
method and the result of Eq. (33) (a) or Eqs. (28), (30) (b), 
respectively. There appears to be only one plot in each figure 
here because of the excellent agreement between these two 
methods. 

used a nested band n = 0.0 with M ~ 0.4 and A = O.lt. 
In both graphs we also plot the curves given by the an- 
alytical results in Eq. (33) and Eqs. (28), (30) respec- 
tively. The (almost) complete overlap of the curves reveal 
that the transfer matrix method captures the low-energy 
bound states, and that these have the same dispersion 
as discussed in the previous sections. Small deviations 
are seen near the Brillouin zone edges. This is expected 
since the Fermi velocity vp vanishes there and the qua- 
siclassical approximation used to derive Eqs. (33) and 
(28) becomes less reliable. 

For the AF/dSC interfaces we find similar agreement 
between the analytical results and the transfer matrix 
method. 



V. THE AF/N/AF JUNCTION 

Discrete quasiparticle bound states below the AF gap 
induced by Q reflection processes at AF/N interfaces, 



exist in a planar AF/N/AF junction analogous to An- 
dreev bound states in SC/N/SC systems^^. Here, we 
present the analytical solution of the problem for the 
simple case As,m <^ vp,x/a, when the sublattice qua- 
siclassical approximation applies well to both the super- 
conductor and the antiferromagnet. Analogous to the 
Andreev equations (20), we formulate Schrodinger equa- 
tions for electrons with the quasiclassical approximation 
in the (kp^kp + Q) representation: 

(^-ivF,x{kp)-^ - II- Ua.kAx) 

-)- am±UcyMF+Q{x) = , 
(^ivp^j;{kp)-^ ~ fi- u„^kF+Q{x) 

+ am±Ucr,kpix) = . (42) 

The relation vp^x{kp + Q) — —vp^x{kp) is taken into 
account in Eq.(42). 

Let be a misorientation angle between two magnetic 
order parameters m± in the right x > d/2 and the left 
X < —d/2 antiferromagnetic halfspaces. Consider the 
magnetizations lying in the yz plane, which is perpen- 
dicular to the interface normal. The quasiparticle spin 
will be not a good quantum number in the case ^ 0, tt, 
so that we have to explicitly introduce spin coordinates. 
We choose the global quantization axis along the magne- 
tization in the left antiferromagnet. Since Eqs. (42) are 
written for the quantization axis taken along the order 
parameters m± > 0, the solutions in the right antiferro- 
magnet should be rotated by the angle 9 around the x 
axis in spin space before matching them with the corre- 
sponding solution in the normal metal region at a; = d/2. 
After the rotation, the low-energy solutions of Eqs. (42) 
in the right antiferromagnetic halfspace can be written 
as 



D 



( cos(0/2)e*^gn("^^-)^+ \ 
isin(0/2)e*«gn(^^^>-)'^+ 



V 



cos(6'/2) 
isin(6l/2) 



D 



/ -isui{0/2)e^^^<'"P--)'^+ \ 
-cos(0/2)e*sgn(i'i=^,r.)¥'+ 



V 



isin(6'/2) 
cos(6'/2) 

I 



cxp 



\Vp. 



(43) 



Here, we have introduced the notation 

^iip±{E) ^J_(^E + n + i^ml - (S + /i)2^ . (44) 

The upper two lines in Eq.(43) describe spin up and 
down amplitudes with the momentum kp and the lower 



two lines show the relative values of spin up and down 
amplitudes with the momentum kp + Q. The solutions 
for the left antiferromagnetic halfspace can be obtained 
from Eq.(43) after the substitutions — > 0, m-^ — s- m_. 

In the normal-metal region —d/2 < x < d/2 the mag- 
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netization vanishes and the solutions of Eqs.(42) take the form 
I 



V TJ-iMF+Qix) 













V 





















Jx{E + tl)/vF.x 



r 




-ix{E + fJ-)/vF,x 



(45) 



In the absence of potential barriers, the quasiclassical so- 
lutions are continuous across the interfaces. Matching 
the solutions (43), (45) at x = ±(i/2, we obtain the fol- 
lowing equation for bound state energies 



\VF,x\ 

2d 



^+{E^) + ip.iEi)±9 + 27:n]-fi , (46) 



where n = 0, ±1, ±2, .... 

The dependence of the bound state energies on the 
particular values m± disappears for low-energy states in 
the almost half-filled lattice. Indeed, under the condition 
\En + n\ <^ m± one can take Lp± ^ tt /2 and find from 
Eq.(46) the following low-energy equidistant spectrum 



E^ikp^ 



K'f,a;(fcf )| 

2d 







2tt -\ 









-M- (47) 



As seen from Eq.(47), i?„+i — i?„ = Tr\vF,x\/d and 
E^{'K)-En{Q) = ±'!T\vF,x\/'^d. The condition \En + ^I,\ < 
m± is valid for sufficiently large thickness of the nor- 
mal metal region d 3> {hvF,x/'ni) ^ ^m, when there can 
be many levels described by Eq.(47). The zero-energy 
bound state appears, in particular, in the half-filled lat- 
tice for = TT when the relative phase of the antiferro- 
magnetic ordering differs by tt in the left and the right 
halfspaces. 

The spectrum (47) qualitatively differs from that in a 
conventional 'particle in a box', i.e. a system of almost 
free quasiparticles confined by impenetrable walls consti- 
tuting a I/N/I junction. The reason for the difference 
is associated with strong correlations between electrons 
with momenta Uf and kp + Q, induced in the normal 
metal region by the antiferromagnets, where equations 
for electrons with kp and kp + Q arc coupled (see Eqs. 
(42)). 

After the substitution m± A the bound state ener- 
gies (46) coincide with those obtained many years ago for 
SC/N/SC systems (at /i = 0) with the phase difference 
6^^ . The reason for this is seen at = where the quasi- 
classical equations Eqs. (42) for electrons with momenta 
kp and kp -\- Q coincide with the Andrecv equations 
with /i = and real A after the substitutions am A, 
'Ua-,kF+Q{x) Va_kp(x). It is important for this prop- 
erty that electrons with the momentum kp + Q possess 
a reverse velocity (similar to holes with the momentum 
kp) compared with electrons with the momentum kp. 



In the presence of a misorientation angle, one can trans- 
form Eqs. (42) into the corresponding Andreev equations 
at /i = with complex A, if one uses a gauge transfor- 
mation after rotating all the quantities in Eq. (42) over 
the angle in spin space. 

In the following, we study the AF/N/AF junction from 
the solution of the BdG equations (4). For the coupling 
constants we have the following simple spatial depen- 
dence 



for 



d 
2' 



V, = for all i, 

and U ^ within the normal region |i| < |. 

(b) 



(48) 
(49) 




FIG. 19: Bound state eigenbands for a 100 AF/N/AF junc- 
tion where the length d of the normal region N is either d = 4 
(a) or d = 48 (b) lattice sites. Parameters: jj, = 0, U = 2.7t 
(a) and U = 2.0t (b). 

In Fig. 19 we plot the eigenbands for the (100) 
AF/N/AF junction when ^ = Q, {/ = 2.7t, d = 4 (a) 
and U = 2.0t, d = A8 (b). In the limit of a long nor- 
mal region, we see that outside the parabola-shaped re- 
gion centered at ky = 0.0, the bands roughly approach a 
set of equidistant sine-shaped bands, in agreement with 
Eq.(47) at 6* = for the (100) orientation. In this hmit 
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FIG. 20: (Color online) LDOS for the 100 AF/N/AF (a) and 
AF/N/ttAF (b) junction with d = i, /j. = 0, and U = 2.7t. 
The four center LDOS scans are in the N region. 



the low-energy LDOS in the N region (not shown) dis- 
plays the expected equally spaced peaks. Note, that in- 
side the parabola region we observe a qualitatively new 
dispersive behavior of the eigenenergies with maximum 
at ky = 0. We expect these deviations from the simple 
result of the quasiclassical treatment to be caused by the 
assumption m <C t used to obtain Eq.(47). In addition, 
the details of the spectrum near points where vp^^ ap- 
proaches zero can always differ from the quasiclassical 
result. 

It is clear from Fig. 19a that we can have very strong 
deviations from the quasiclassical result when d ^ a and 
m ^ t. In particular, within the quasiclassical framework 
the dispersion of the bound state energies is entirely as- 
sociated with the momentum dependence of the Fermi 
velocity VF,x{kF)- In contrast, Fig. 19a displays a more 
complicated momentum dependence of the energy lev- 
els with numerous extrema and additional peaks in the 
associated LDOS. 

We have also studied a spin-pi junction, a AF/N/ttAF 
junction with 9 = ir, where the magnetization has gained 
an extra tt phase shift when crossing the normal region 
N. This is similar to the 7r-phase shift in the stripe phase 
of the cuprate superconductors^^"^^. For the AF/N/ttAF 
junction the subgap bands are shifted compared to the 
AF/N/AF configuration. As a consequence, the associ- 
ated LDOS will reveal whether the junction exhibits the 
TT phase shift or not. This can be seen in Fig. 20 which 
shows the LDOS associated with Fig. 19a and the corre- 
sponding AF/N/ttAF junction in Fig. 20b. The surpris- 
ingly large number of peaks seen in Fig. 20 agree with 
the strongly dispersive bands as can be verified from the 
stationary points in Fig. 19a. The presence of low-energy 
states seen in the AF/N/ttAF junction in Fig. 20b cor- 
relates with predictions of the quasiclassical approach. 
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FIG. 21: (Color online) Eigenbands and the associated LDOS 
for the 110 AF/N/AF (left column) and 110 AF/N/ttAF 
(right column) junction with d = 8, /i = 0, and U = 2.7t. 
Note the low-energy state generated in the spin-Tr junction. 
In the LDOS plots, the middle eight scans are in the N region 
whereas the upper and lower three shows the LDOS upon 
moving into the AF. 

which, however, does not apply directly to the case of 
very thin normal metal region d = A. 

Fig. 21 shows the bands and LDOS for the 110 
AF/N/AF (left column) and AF/N/ttAF (right column) 
junction, respectively. In accordance with Eqs. (46), 
(47) , the orientational dependence of the spectrum is as- 
sociated with the Fermi velocity dispersion. Indeed, in 
the limit of large d, the bound state spectrum consists 
again of equidistant states like in the (100) case. For 
the 110 geometry we also find that the main qualitative 
difference between the AF/N/AF and AF/N/ttAF junc- 
tions is associated with the presence of almost zero energy 
states in the latter case. 



VI. CONCLUSIONS 

We have performed a theoretical study of interfaces 
and junctions involving antiferromagnets and supercon- 
ductors or normal metals. This was presented both in 
the framework of quasiclassics and self-consistent numeri- 
cal solutions of the relevant Bogoliubov-de Gennes equa- 
tions. Where comparison is appropriate, we found full 
agreement between the two methods. In particular, we 
investigated for formation of bound states near (100) and 
(110) interfaces between antiferromagnets and s- or d- 
wave superconductors. We calculated their dispersion, 
their influence on the proximity effect, and the associated 
modifications of the LDOS near the interface regions. In 
addition we discussed the crossover between Q-reflection 
and conventional specular reflection as a function of the 
potential barrier at the interface. In future work we plan 
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to investigate the role of the bound states on the Joseph- 
son current in SC/AF/SC junctions where the low-energy 
bound states can be expected to generate e. g. unusual 
temperature dependence of the critical current. 
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